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ABSTRACT 

Context. Trails of dust inside galactic bars are easily observable in external galaxies. However, information on the dust lanes of the 

Milky Way bar is harder to obtain due to our position within the Galactic disc. 

Aims. By comparing the distribution of dust and gas in the central regions of the Galaxy, we aim to obtain new insights into the 

properties of the offset dust lanes leading the bar's major axis in the Milky Way. 

Methods. On the one hand, the molecular emission of the dust lanes is extracted from the observed CO l-b-V distribution according 

to the interpretation of a dynamical model. On the other hand, a three dimensional extinction map of the Galactic central region 

constructed from near-infrared observations is used as a tracer of the dust itself and clearly reveals dust lanes in its face-on projection. 

Comparison of the position of both independent detections of the dust lanes is performed in the (/, b) plane. 

Results. These two completely independent methods are used to provide a coherent picture of the dust lanes in the Milky Way bar In 

both the gas and dust distributions, the dust lanes are found to be out of the Galactic plane, appearing at negative latitudes for / > 0° 

and at positive latitudes for / < 0°. However, even though there is substantial overlap between the two components, they are offset 

from one another with the dust appearing to lie closer to the b = 0° plane. 

Conclusions. Two scenarios are proposed to explain the observed offset. The first involves grain destruction by the bar shock and 

reformation downstream. Due to the decrease in velocity caused by the shock, this occurs at lower z. The second assumes that the gas 

and dust remain on a common tilted plane, but that the molecular gas decouples from the Milky Way's magnetic field, itself strong 

enough to resist the shear of the bar's shock. The diffuse gas and dust remain coupled to the field and are carried further downstream. 

This second scenario has recently been suggested in order to explain observations of the barred galaxy NGC 1097. 
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. 1. Introduction In this letter, we combine the results from gas dynamics and 

I — ^ ,,. ,.. ..TT- ». ,-. Stellar reddening studies in order to provide a descr iption of the 

r- , Barred galaxies are ubiquitous in the Universe. Many of these ^^^^ j^^^^^ j^^^. ^^ ^j^^ ^^^,^ ^^.^^ ^^.^ P^- r^^^ interpreted 

Vi : are observed to have dust lanes - trails of gas and dust, seen in ^j^^ j^j ^^^ ^q longitude-velocity diagrams of the Galactic disc 

> extinction, leading the bar s major axis. Evidence for the bar na- ^^. ^elf-consistent N-body and smooth particle hydrodynam- 

• -H ^ t ure of the Milky Way s bulge has been around for over 40 years .^^ simulations of the stellar and gaseous components. Later, 

^ , d^de VaucouleursI [19^1). However our position in our Galaxy Marshall iH] (gOOi) used a galactic stellar population synthe- 

^ , does not give us a clear view of the inner regions. Photometric ^.^ ^^^^j ^^ ^^^^^^^ ^j^^ ^j^^^^ dimensional extinction informa- 

■ - - studies are hampered by the high source confusion and severe in- ^^^^ ^^^^ ^^^ .^^^^^.^^ ^^^U^^. observations. Both studies were 

terstellar extinction whereas gas kinematic studies are rendered ^^le to identify the Milky Way's dust lanes. In the following, we 

difficult due to line of sight crowding of the gas emission lines. ^ ^^^ ^^^^ ^^^ ^ ^^^ ^^ distributions of the 

Nevertheless many studies have obtained information on ^^^^ ^^^^^ ^^^^^^. ^^^^ ^^^^^ ^^^^.^^ 
the Milky Way s dust lanes using a variety of techniques. 
ICalbet et aTl d 19961) detected a longitude assymetry at 10° < 

|/| < 20° in the DIRBE s urface brightness maps which they at- 2. Dust lanes from 3D extinction maps 

tributed to the dust lanes. iBabusiaux & Gilmora (l2005h used red , , j , 

clump stars to determine the distance as a function of longitude [Marshall et al. | (12006D created a three dimensional map of the 

to the bar and using their observations concluded that there is 1""^"" G^l^'^y ^^i^g *e Two Micron All Sky Survey(2MASS) 

no near-infrared extinction due to the dust lanes in our Galaxy. Po^"! source catalogue (PSC) ( |Skrutskie et al. | |2006|) and the 

Other studies have based their interpretation of the observa- Stellar Population Synthesis Model, developed in Besangon 

tions on models in order to s tudy the chemistry in the G alactic <Robin et al. 2003, which we will simply call "Galactic model"). 

centre (GC) and dust lanes (iRodriguez-Fernandez et al.ll2006t For a full descrJEtion of the 3D extinction method, please see 

iNagai et alJbOOTl) as well as the physical processes underway [ Marshall et al. |( |2006|). 
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in the dust lane shocks CTszt 2006). However, the results may Contrary to [Marshall et al. | (12006!), the method used in this 



depend on the assumed geometry for the dust lanes as no coher- ^tudy does not have a fixed angular resolution but one that varies 

ent picture currendy exists for them ^^ ^ function of the stellar density. The extinction information 

for the GC region presented here has an angular resolution that 

* Free-time astronomer varies from 3.75' (high density) to 15' (low density). 
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Fig. 1. Application of the extinction method to the line of sight 
I - 356°, b = 0.75°. The crosses represent the total extinction 
as a function of distance (left hand scale) and the solid red line 
is the corresponding differential extinction along the same line 
of sight (right hand scale). The differential extinction is assumed 
proportional to the dust density. 



The G alactic model used here differs slightly from that de- 
scribed in lRobin et al.l (l2003l) . as a result of more recent adjust- 
ments using 2MASS data. Here, the stella r density of the bulg e 
is modelled using the sech function from lFreudenreichI (Il998h . 
It is boxy and prolate, with axis ratios 1 : 0.3 : 0.25, and it is 
assumed to form an angle of 20° with the Sun - GC direction, in 
line with recent estimates (lGerhardll2002h . 

The three dimensional extinction distribution for an example 
line of sight is shown in Figure [1] Each cross represents an es- 
timation of the total extinction at the given distance. In order to 
visualise the three dimensional dust distribution we will present 
a map of the differential extinction, which is the derivative of the 
extinction with respect to distance. This will be proportional to 
the density of the big dust grains (15n m < radius < llO nm) as 
they dominate the extinction at 2 yum jDesert et al.lll990 ). The 
differential extinction for the example line of sight is shown as 
the solid red line in Figure[T] Under ~ 2 kpc and at very large dis- 
tances, where the distance-extinction points are spaced far apart, 
it is difficult to identify any underlying structure. However, the 
technique is well suited to studies of the Galactic bar : the ex- 
tinction due to the dust lanes for this line of sight can be seen at 
r ~9.5 kpc. 

The distribution of dust in the Galactic plane, towards the 
GC, is displayed in Figure|2l The elongated structure seen in ex- 
tinction, running along the stellar bar (shown by the solid line 
in Fig|2]) but with a slightly different angle, is due to the dust 
lanes. They can be seen to precede the bar in its clockwise mo- 
tion (as seen from the NGP), in agreement with observations of 
dust lanes in external galaxies. The near-side dust lane becomes 
intertwined with the molecular ring near x = 6 kpc, y - 1 kpc 
and the far-side dust lane reaches out to distances where the ex- 
tinction map becomes unreliable or non-existent. 

3. Dust lanes from molecular gas kinematics 

Figure [3^ shows the observed / - V distribution of '^CO in the 
Galactic bar region. The black contours en close the two features 
which have been identified bv lFuxl (1 19991) as the gaseous signa- 
ture of the dust lanes along the Galactic bar. The I > 0° feature, 
also called "Connecting arm", is the trace of the dust lane in the 
near-side part of the bar, and the I < 0° feature is the trace of the 
far-side dust lane. 

Because of the distance indetermination inherent to l-b-V 
observations, one cannot entirely isolate the emission of these 




Fig. 2. The average density of absorbing matter in the Galactic 
bar region for |z| < 300 pc, as seen from the NGP. The quantity 
plotted is the differential extinction (in units of magnitudes per 
kpc in the K, band). The GC is indicated by the square sym- 
bol, the solid line shows the orientation of the modelled stellar 
bar, and the ellipse shows the region selected, by inspection, to 
contain the extinction arising from the dust lanes. The Sun is at 
(0,0), outside the figure. The dashed lines indicate galactic lon- 
gitudes I = -5° and / - 10° which approximately delimit the 
ends of the dust lanes as seen in CO. The white band indicates 
where the extinction estimates do not fully sample the bar due 
to the low completeness (high source confusion) of the 2MASS 
PSC. 



dust lanes from foreground and background emission. The fol- 
lowing describes our approach to select most of the dust lanes 
CO emission while excluding most of the unrelated emission. 
The first of the adopted selection criteria restricts the emission 
to the I- V regions highlighted in Figure[3^. This criterion avoids 
in particular the emis sion of the |/| ^ 1.5° nuclear ring/disc and 
of Bania's Clump 1 feanial 1977b at (/, V) ~ (-5°, 100 kms"'). It 
also avoids the emission of the velocity elongated features near 
/ = 3° (Bania's Clump 2) and / - 5.5°, which may arise from 
molecular clouds being absorbed by the near-side dust lane dFuxl 
1999; Liszt 2006; Rodriguez-Fernandez et al. 2006). 

The second criterion concerns the gas so far selected in the 
near-side lane / - V region. The / - b distribution of this gas 
(Fig. [3}5) reveals independent emission within 3.125° < / < 
10.5° and -0.125° < b < 0.625° (black box in the figure) that 
forms two distinct I - V features (Fig. [3]:). The nearly constant 
velocity orientation of these features contrasts with the global 
/ - V inclination of the near-side lane and thus suggests that they 
do not belong to this lane. The gas within the box in Fig. [3}5 
is therefore discarded. Finally, the third selection criterion con- 
cerns the gas selected in the far-side lane I - V region. The b-V 
distribution of this gas (Fig. [3}l) separates fairly well the emis- 
sion of the dust lane from the bright emission of external spiral 
arms, near V = Okms ' and centred at b a 0°, and of the 3- 
kpc arm, at (b, V) ^ (0°, -75 km s"'). Hence we also exclude all 
emission between the two lines drawn in Fig. [3}l. It should be 
noted that the remaining emission at / ^ 8.5° has characteristics 
similar to those of the features seen in Fig. [3}; and thus may also 
not originate from the dust lanes. 
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Fig. 3. Selection procedure of the '^CO e mission assoc i ated w ith 
the Galactic bar dust lanes, based on the iDame et aTl ( 1200 ll) l- 
b-V data cube, a) Longitude-velocity diagram of the CO gas 
within \b\ < 2.5°, with the black contours indicating the near-side 
(Connecting arm) and far-side lane regions set by our first selec- 
tion criterion, b) Longitude-latitude diagram of the gas within 
the near-side lane region, c) Longitude-velocity diagram of the 
gas enclosed by the black box in the former l-b plot; this con- 
tribution is rejected according to our second selection criterion. 
d) Longitude-latitude diagram of the gas in the far-side lane l-V 
region; our last criterion excludes the emission between the two 
displayed lines. 



4. Results and discussion 

The sky distribution of the extinction due to the dust lanes is 
displayed on the top of Figure H] The quantity plotted is the in- 
tegrated extinction owing to all the matter inside the ellipse of 
Fig. 121 We choose sky projections to compare the extinction of 
the dust lanes with their CO emission to avoid any modelling 
required to deproject the raw CO data. In the bottom of Figure|4] 
is the I - b distribution of the gas as selected in Secl3] The tilt 
of the dust lanes relative to the Galactic plane is obvious in both 
the gas and the dust components, with the negative longitude 
(far-side) part sitting above the plane and the positive longitude 
(near-side) part below the plane for both components. Also, both 
components are well contained w ithin |/| < 10° , show ing that we 
are tracing a different structure to lCalbetet alj ( ll996l) . 

The mean plane of the dust/gas as a function of galactic lon- 
gitude is given by b{l) ^'Zbb- X{1, b)/Y,h X(l, b) where X{1, b) is 
either the extinction due to dust or the intensity of the CO emis- 
sion emanating from the dust lanes. The calculation of the mean 
plane for the dust uses lines of sight where Aks > 0. 15 (indicated 
by the solid contour in the top image of FigHJi, in order to avoid 
low density lines of sight unrelated to the dust lanes. The mean 
planes of the dust and the gas, calculated this way, are shown 
in Figure |5] Both the gas and the dust follow the same trend, 
however the dust is seen to be displaced from the gas. 

The magnitude of the apparent displacement is larger at pos- 
itive longitudes than at negative longitudes. The removal from 
the CO observations of the zone highlighted in FiglSj^ may bias 
the mean plane to more negative latitudes but, as the bulk of the 




Fig. 4. Longitude-latitude distributions of the near-infrared ab- 
sorbing dust lanes detected in the Galactic bar by Marshall et al] 
( l2006h ( top), and of the co rresponding '^CO emission derived 
from the lDame et alJ ( 200lh d ata assuming the l-V feature inter- 
pretation given bv lFuxl(ll999l) (bottom). Both images are filled 
contour plots with 7.5'resolution. The Galactic b - Q° plane is 
indicated by the dashed line. The CO emission from the central 
molecular disc has not been included (see Sec|3j, and the dust 
map excludes the area masked out in Fig|2] The contour rep- 
resents Aks - 0.15 mag. Both the gas and dust are seen to be 
tilted with respect to the Galactic plane, with the negative longi- 
tude (far) side sitting above the plane and the positive longitude 
(near) side sitting below. 



dust lane CO emission arises from outside of this zone, the ef- 
fect would be small. Likewise, the contribution of the dust com- 
ponent outside the dust lanes has been minimised using the high 
extinction lines of sight chosen using the ellipse of Sec|2] In fact, 
the difference in apparent displacement is largely an effect of 
perspective, as the near side of the bar lies at positive longitudes. 

These results suggest that, even though the large scale mor- 
phology of the dust lanes is similar for the two components stud- 
ied here, there is an additional effect which causes their /, b pro- 
jections to be offset from one another This could be due to an 
offset in z, a separation in the x, y coordinates of a tilted plane or 
a combination of the two. 

Offsets have been detected between dust and gas in external 
galaxies seen face on, where azimuthal offsets are easier to dis- 
cern than offsets in z- For example, in M83 the dust extinction 
and CO em ission in the eastern spiral arm are offset by as much 
as 700 pc (lLord&Kennevlll99U iRandet all [19991) . Three ex- 
planations are offered to explain the offset : 1) The CO is heated 
by UV radiation from young stars; 2) The CO is heated by low- 
energy cosmic rays; 3) The components react differently to the 
spiral density wave. In the case of the Milky Way bar, the first 
option is unlikely as star for mation along the bar is likely harn - 
pered by the high shear there dRodriguez-Fernandez et al.l2006l) . 
and stars for med at the ends of the ba r are not expected to mi- 
grate along it dCole & W einberg 2002). Cosmic ray heating may 
be responsi ble for increased gas tem peratures in the centre of 
the Galaxy dYusef-Zadeh et al.ll2007h but these tem peratures can 
also be explained by shock waves dLis et alj|2001h . In the third 
case, the diffuse gas is trapped in a thin dust lane at the shock 
front and the molecular clouds form a broad ridge on the spi- 
ral arm. However this does not seem to correspond to our ob- 
servations, where it is the dust that is seen to be more broadly 
distributed than the molecular gas (Fig|4]i. 

Recently, Bec k et al] y005) observed the barred galaxy 
NGC 1097 in total and polarised intensity at A =3.5 cm and 
A =6.2 cm. They found that around the bar the total intensity in- 
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creases along the shearing shock, however the increase of the po- 
larised emission occurs further downstream and is aligned with 
the bar's dust lanes. They suggest that the large scale magnetic 
field is only weakly affected by the shock and that the molecular 
clouds become decoupled from the magnetic field and are thus 
swept up by the bar's shearing shock. The dust, coupled to the 
magnetic field, thus finds its way ~ 400 pc downstream from 
the CO. Furthermore, they suggest that this may occur in most 
disc galaxies where large scale shocks trigger the formation or 
collapse of gas clouds. Our observations are consistent with this 
scenario, assuming that the gas and dust of the dust lane s orbit 
on a common and tilted plane. However, B eck et al] ('2005) were 
not able to obtain direct observations of the molecular gas - dust 
separation in this galaxy. In order to detect an offset of a few hun- 
dred parsecs between CO and dust in an external galaxy would 
require high resolution and high sensitivity CO observations of 
a nearby barred galaxy with prominent dust lanes. 

Instead of considering an offset on a common plane, it is pos- 
sible that the observed offset is in effect a separation in z- As the 
gas and dust hit the shearing shock, it is possible that a fraction 
of this gas and dust fall towards the Galac tic plane. Two possible 
mechanisms could be considered. Firstly. lOomez & Coxl (|2004|) 
show using 3D MHD simulations that gas hitting a spiral arm 
shock can be deflected over the arm. Although there has been 
no study of this kind done on the bar shocks, there is also no 
evidence indicating that the gas and dust should be confined to 
the same plane. Secondly, the molecular gas and dust enter the 
shock front and merge with the dust lane at a given height above 
the Galactic plane. Due to the gravitational potential, the post- 
shock material eventually falls onto the Galactic plane creating 
a vertical distribution of azimuthally comoving gas and dust in 
the dust lane. 

These mechanisms could explain a vertical distribution of 
gas and dust at the shock fronts but not a separation between 
them. However, at the shock front the big grains of the dust 
(responsible for the NIR extinction) are shattered jJones et alj 
11994) . creating many smaller grains (responsible for UV extinc- 
tion). As the gas moves away from the shock front, the increased 
density speeds up the rate of coagulation between grains and, in 
particular, small grains on large ones. If some of the gas and dust 
are heading towards the Galactic plane, the lack of CO at low z 
could be explained by two processes: 1 ) accretion of gaseous CO 
on to grains a s an ice mantle, which is efficient for Aks ~ 0.35 
(lDrainell2003l) : 2) dissociation of CO due to the increase in UV 
photon flux resulting from the decrease in the small grain abun- 
dance. The dust layer would be found at systematically lower z 
than the CO, as can be seen in Figs.|4]and|5] 

iFuxl (1 19991) represented the tilts of the CO and HI by a 
straight line with an angle of 6 ~ 4.5° at their azimuth of 25°, 
while mentioning that this representation is an oversimplifica- 
tion. The tilt of the dust component is not as easily represented 
by a straight line (Fig|5]). Its value is highly dependent on the 
segment of the dust lane which is fitted and the interpretation 
of its value depends on the scenario used to explain the offset 
discussed above. 

Finally, there are uncertainties in the dust position, namely 
in its distance along the line of sight. This uncertainty is es- 
timated to be about 500 pc in the bar (Fig. [1]). Moreover the 
computation of the 3D extinction map made use of the Galactic 
model, hence it depends slightly on the parameters of this model 
(Marshall et al. 2006). This dependence might cause a system- 
atic error of a few 100 pc in the distance of the dust along the 
line of sight, but not on its (I, b) position. It is expected that this 
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Fig. 5. Position of the mean plane of dust (crosses) and molecu- 
lar gas (diamonds) along the dust lanes in the Galactic bar The 
b - 0° plane is shown by the dashed line. The same large scale 
tilt is observed for both the gas and dust, however the dust ap- 
pears closer to the b -Q° plane. 



error would affect the positive and negative longitudes similarly, 
having no influence on the tilt of the dust lanes in the (/, b) plane. 



5. Conclusions 

We have identified the dominant dust lanes of the Galactic bar 
both by their near-infrared extinction and by their CO emission. 
The same large scale behaviour is observed in both components, 
that is a tilt going from /? < 0° at positive / to /? > 0° at negative 
/. In the /, b plane projections, the two components are seen to be 
displaced from one another This offset may be due to dust de- 
struction at the shock front followed by reformation and/or ice 
mantle accretion downstream and at lower z, or the molecular 
gas becoming decoupled from the magnetic field which resists 
the shearing shock of the Milky Way's bar and carries diffuse 
gas and dust further downstream on a tilted plane. Further obser- 
vations of our Galaxy and external galaxies, together with the- 
oretical work are necessary in order to unambiguously interpret 
these observations. 

Finally, the coincidence between the two independently re- 
covered gas and d ust d i stribu tion along the bar reinforces the 
interpretation from lFuxl d 19991) regarding the I -V traces of the 
dust lanes in the observed HI and CO kinematics, and it adds 
credibility to the Galactic model and the use of it to determine 
the three dimensional extinction distribution. 
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